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Abstract 

The effect of the angular momentum density of a gravitational 
source on the times of flight of light rays in an interferometer is an- 
alyzed. The calculation is made imagining that the interferometer is 
at the equator of an axisymmetric steadily rotating gravity source. 
In order to evaluate the size of the effect in the case of the Earth a 
weak field approximation for the metric elements is introduced. For 
laboratory scales and non-geodesic paths the correction due to the 
angular momentum turns out to be comparable with the sensitivity 
expected in gravitational waves interferometric detectors, whereas it 
drops under the threshold of detectability when using free (geodesic) 
light rays. 
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1 Introduction 

The famous Michelson-Morley experiment does not require any explanation 
regarding its nature and the crucial role that history reserved to it is well 
known. It has been analyzed in any respect in the early days of relativity and 
discussed also on its fundamental meaning Since then it has been assumed 
that no anisotropy can be revealed until the frontier of special relativity is 
crossed. Only in a few cases anisotropics deriving from general relativistic 
corrections were considered [0]0, but only caused by the gravitational red 
shift in non-horizontal arms of the interferometer; Schwarzschild-like correc- 
tions do not produce any effect in the horizontal plane. 

However, if the source of the field is rotating as it is the case for the 
Earth, the situation in principle changes. This means that a tiny anisotropy 
can legitimately be expected, depending on the angular momentum of the 
source. 

On the other hand the search for measurable effects of the angular mo- 
mentum of the gravitational field is always active in order to add a new 
direct verification of the consequences of general relativity. The only positive 
result at the moment concerns the precession of the nodes of the orbit of 
the LAGEOS satellite [Q (Lense-Thirring effect [^). In the next few years 
the space mission Gravity Probe B (GPB) is planned to fly carrying gyro- 
scopes which should in turn verify the Lense-Thirring effect too finally a 
series of different possibilities connected both with the Sagnac effect and the 
gravitomagnetic clock effect have been considered 

The present paper will fix the general formalism to verify a possible influ- 
ence of the angular momentum density of the Earth on a Michelson-Morley 
type experiment. Numerical estimates will show that the effect is quite small 
in any case, however using non-geodesic light paths it could turn out to be 
comparable with what people are expecting and planning to measure with big 
size interferometric gravitational waves detectors like LIGO and VIRGO 



2 Preliminaries 

The Michelson-Morley experiment is an interferometric measure and uses 
light, let us then start from a generic null line element in polar coordinates 
and within an axially symmetric static field originated by a central body 
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endowed with an angular velocity Q = d(f)/dt: 



= gudt" + 2(7, 



grrdr"^ + geedO"^ + 



(1) 



The (yf's are of course the elements of the metric and are independent both 
from time t and from 0. In weak-field approximation the explicit form of the 



metric can be 11 
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where we introduced the parameters a = J/Mc (J is the angular momentum 
of the source, M is its mass and c is the speed of light), fi = 2GM/c^ 
(Schwarzschild radius of the source). Now let us consider r = constant world 
lines only. This choice corresponds to limiting the study to light beams 
contained locally in a "horizontal" plane (actually this would require a wave 
guide locally shaped as a constant gravitational potential surface). The null 
(non-geodesic) world line becomes: 



= gudt" + 2g, 



(3) 



If we suppose to place our interferometer at the equator {6 = n/2) and 
provided its arms are not too long, the metric (§) (first order in a/ r, fi/r and 
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Consequently for short enough excursions in the "horizontal plane" we can 
assume, at the lowest order in 6*, that, for light, (f) and 9 variations are 
approximately proportional to each other, so: 

\de\ = x\d<p\ (5) 

where x is a constant. 

Suppose now that the interferometer arms are stretched one in the North- 
South direction and the other in the East- West direction. Taking into account 
the fact that the Earth reference frame, where the interferometer is at rest, 
is indeed rotating, the coefficient x 'wiH depend on the angular speed of the 
Earth. 

Now solving (^) for dt we obtain 

-9t4> ± J 9t^ - 9tt9eeX^ - 9tt9^<p , , 

dt = d(p (6) 

9tt 

Supposing our light beam starts from a point on the equator at = 
and moves northward, it will be 

= l-X<l> (7) 
Of course in the case of an East- West beam it is x = 0. 



3 Times of flight of non-geodesic light beams 

In proximity of the equator the first factor in the right hand side of does 
not depend on 6; thus 

-9t<j, + J 9^4 - 9tt9eeX^ - 9tt9<l><t> , . . 

tN = 91 (8) 

9tt 

t^ is the time of fiight to reach the northern mirror and (pi is the angular 
coordinate of the event; the drift of the beam is naturally in the prograde 
sense. 
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The world line of the mirror (initially at 6' = | — $ and (f) = 0; here $ 
represents the angular stretch of the interferometer arm) is: 

tN = h/^ (9) 

Q is of course the angular speed of the Earth. 

(H) and (I) allow to deduce an expression for x- 



1 / gtt + 2gt^n + g^^n^ 
^=0V ^0 

Actually at the lowest order in 6 (|Ty) does not contain 6 itself and conse- 
quently on an r = i? = constant surface x is a constant. 

From d^) we see that to span the South-North angular distance one 
travels eastward by the angle 

^i = - (11) 

X 

Returning to (^) and using ([T^): 



IN 



gtt + Sfiff^fi + g^^^l'^ 



Considering the North-South way back to the source we see from (|^) that 
it is 

— fl't0 + \j9t4,~ QttQeeXs ~ QuQcj,^ 
ts = 02 (12) 

gtt 

where 02 is the Earth's rotation angle between the reflection and the arrival 
back at the source. It must be 

ts = 02/1^ (13) 

([T^ ) and (pT^) give xs = X iii (FOD; then it is ts = tN, 02 = 0i. Finally the 



total time of flight South-North-South is 



tsNS =tN+ts = 2 

V gtt + 2gt^il - 
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To proceed further we recall the explicit expressions for the g^s, given in 
d^): now, posing ^ = l/R where / is the length of the arm of the interferom- 
eter and R is the radius of the Earth, one has approximately: 

All further corrections in 6, i.e. $, are indeed quadratic and multiply the 
other small terms, thus resulting much smaller than them. 

The next step is to consider the time of flight along the East- West arm 
of the interferometer. 

From formula with 6* = 7r/2 = const we have 

-9t<p + J 9t^ - 9tt9<p4, 

tE = (14) 

9tt 

The world line of the eastern end mirror, assuming equal length arms, is 
= $ + f2t which means also 

(j)E = '^ + ^tE (15) 

Combining (p!4[) and ([T5|) one has 

tE = ^\ \ <^ (16) 

9tt + 9tS - J [9t4> - 9u9<P4>)^ 

Now for the way back we have 



9tcp + 9t^ - 9tt9H , . 

tw = {(PE - (Pw) (17) 

9tt 

(pw is the angular coordinate of the source at the arrival time of the reflected 
beam. It must also be 

(l)W = ^ (tE + tw) (18) 
(|T8|) and ( [T7| ) together allow for the calculation of tw- 



tw = ^W^^^^$ (19) 

9tt + 9t<p^ H 



The total West-East-West time of flight is 

\J 9t4, - 9ttgH 

twEW = tE + tw = 2 — - — „ $ (20) 

gtt + 2gt^il + iVg^,^ 

The difference in the time of flight along the two arms at the lowest order 
weak field approximation is 

= twEW - tsNS - (21) 



4 Geodesic light beams 

The situation for free, i.e. geodesic, light rays is different from the description 
given in the previous section. Now we start from the remark that in the 
equatorial plane the bending of the light rays is lower than the curvature of 
the circle along which the mirrors of the interferometer move. Actually in 
the zeroth order of approximation the time of flight of light between the two 
end mirrors of an arm of the interferometer is deduced from the length of the 
chord subtended to the appropriate arc of the mirrors circumference. It is 

to = 2- sin ^ (22) 

where (/)o = $ ± fit (-|- in the prograde path, — in the reverse trip). 

Taking into account the effect of the mass M and the angular momentum 
density a, we expect a deviation from the straight line, which can be expressed 
in terms of the space curvature of the light beam k = 1/p {p is the radius of 
curvature) . 

A further approximation can be to use the average curvature of the path 
between the two ends of the interferometer arc. The length of the intercepted 
beam would then be 

(23) 

where ip is the angle subtending the moving interferometer arc, as seen from 
the curvature center. When the space trajectory of the light rays is not 
contained in the equatorial plane, we expect it also no more to be plane at 
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all, however reasonably the non-planarity corrections will be smaller than the 
other corrections we are introducing. 

The chord subtended to the arc ( p3D is of course the same when seen from 
the origin of the reference frame; in the equatorial plane this gives 

p sin — = R sin — 
Now is slightly different from the former 0^: 

= 0o ± 

where 5t is the correction to the time of flight induced by the curvature of 
the trajectory. 

Extracting -0 from (^: 

, o . fR . A 
lb = 2 arcsm — sm — 

and the time of flight becomes 

P , ■ f R ■ <P\ 

t = —w = 2— arcsm — sm — 

c c \p 2) 

More explicitly in terms of 5t and using (p2D: 

e, , , .P . (R . 0o±^Wt\ R . <t>o 

ot = t — tn = 2— arcsm — sm — 2 — sm — 

c \p 2 ) c 2 

Reasonably it is p >> i? and all angles are small (< 10"^ rad, which is the 
angle subtended under aim arm, from the center of the Earth). This allows 
for series developments up to the lowest meaningful orders, thus 

Solving for 5t 

1 

5t=——- (24) 
24 cp2 ^ ' 



Now we need an explicit expression for p. A standard approach |Tl 
moves from considering the right hand side of divided by (where A is 
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an affine parameter) as the Lagrangian of the hght ray. The cychcity of the 
t and coordinates leads to the constants of the motion 

dt d(f) 
E = 9u^^+9t,^ 

dt d(j) 



wherefrom 



9t<t>E - guL 



9t<p - 9ttg,j,^ 



d\ 

dt -g<t><t>E + gt^L 



d\ 



2 

9t<i> ~ 9ttg4>(i> 



Using the Lagrangian in the equatorial plane one obtains also 

( drV _ 2ELgt^ - E'^g^^ - L'^gu 



A^J grr {gttgcp^ - gt^ 



(25) 



and 



^^,^2 ^ 9ttg<p<i, - gt4, ^jgt^ - g^^ - fgu ^25) 
9rr {jgtt - gt<t>f 



A ' denotes differentiation with respect to and j = L/E. Introducing 
the variable u = 1/r one has 

^^,^2 ^ ^^ gtt{u)g^^{u) - gl^{u) 2jgt^{u) - g^^{u) - fgaju) ^^7) 
9rr{u) Ugttiu) - gt^{u)f 

Now using (|2D (^7]) , in the lowest order in the gravitational parameters, be- 
comes 



iu'f = ^-u^ + iiu^ + 2^fia + ( 4 - 2mM u^a^ (28) 

j2 j6 \j2 



Differentiating with respect to we end up with the differential equation 
u + u = -au H — -aa — Au a -\ — -ua (29) 

2 ji j2 



9 



Now coming to the local curvature of the rays in the equatorial plane, it 
is convenient to use a Cartesian reference frame in that plane posing 



y = \l-g<t>^ cos< 



In terms of u (which is a function of (p along the trajectory) and using the 
appropriate approximation: 



then 



Finally 



y = — \/(l + a^M^) cos< 
u * 

X = — v(l + a^M^) sine; 
u * 



dy u' cos (f) + {1 + u'^a^) u sin (f) 
dx u'sin</) — (1 + n^a^) ucos<^ 

1 dw 



A{(t)) + B (0) a' (30) 
p dx 

where A (0) and E (0) are rather complicated functions of m, m', u" and 0. 

The formula (|5D|) can be explicitly written in a convenient form if the 
curvature is calculated at the point of closest approach of the ray to the 
center of the Earth (maximum value of u which we are calling Um)- There 
we expect m' = and can decide that = 0. Reasonably the average value of 
the curvature along the path of the light differs from the value at the closest 
approach by small corrections. 

The curvature value we shall use is then 



U'L +Um) + - [Um + 3m^) U^O^ 



p ^ " 2 

Actually the last term is 12 orders of magnitude smaller than the first, so in 
practice 

- = - (C + Mm) (31) 

Now be obtained equating (^) to 0. At the lowest approximation 

order it is 

_ 1 

Mm ~ 

3 
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Introducing this result into and (PT| ) produces 




9/i 0/^2 
— 7 — 3— a 

4^-4 ^-5 



Now back to the time difference (^41). A further crude approximation can be 
to put j ~ R: 

1 /3 /i /9 ^a2^ 



(32) 

24 c \^4^ y ^ ^ 

The correction, as we see, is neghgibly small: ~ 10~^^ s for the pure mass 
term and ~ 10"^'' s for the a? term. 

Out of the equatorial plane we cannot reasonably expect the situation 
to be different. Of course the difference in flight times along the two arms 
cannot but be less or at most equal in the order of magnitude to (^). In 
practice this means that for free light rays in the terrestrial environment the 
straight line approximation is fairly adequate and the time of flight difference 
has the typical value of the Michelson-Morley experiment. 

5 Discussion and conclusion 

The result (^) is obtained under the assumption that a physical apparatus 
(bidimensional wave guide) obliges the light rays to move along constant 
radius paths. Were this possible the order of magnitude estimate at the 
surface of the Earth for 1 m long interferometer arms would be: 

At ~ 10~2° s (33) 

This effect is purely rotational and rather small but not entirely negligible. 
Should the interferometer rotate in the horizontal plane, the time of flight 
difference At would alternatively change of sign displaying an oscillating 
behavior which in principle could be detected. In general ([33| ) fixes the scale 
for effects on the Earth. 



A way to increase the value of (|33|) would be to have multiple reflections 



of the light beams happen along each arm of the interferometer before the 
actual interference is measured. This is what happens in Fabry-Perot type 
interferometers. Here both arms should be equipped with such devices, it 
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would then be easy to have the hght rays to be reflected back and forth 
for instance 10^ or 10^ times, before the measurement. Since the effect we 
are looking at is indeed cumulative, this fact would bring the time difference 
between the two paths to 10~^^ — 10~^^ s, which corresponds, for visible light, 
to 10~^ — 10~^ fringe in the interference pattern. This would indeed be a 
huge effect: 10~^ — 10~^ fractional change of the signal intensity. 

It is remarkable that the obtained numeric value compares with the ex- 
pected phase (and time) shifts in the gravitational wave interferometric de- 



tectors now under construction, as LIGO and VIRGO [§])110i- There indeed 
a sensitivity is expected, in measuring displacements, of the order of 10~^^ m 
which corresponds to a time of flight difference 4 orders of magnitude lower 
than (^) and consequent much higher sensitivity. 

Of course our effect as such is a static one, producing a DC signal and it 
would be practically impossible to recognize its presence in the given static 
interference pattern. On the other hand the spectacular sensitivity of gravi- 
tational wave interferometric detectors is obtained at frequencies in the range 
10^ — 10^ Hz. To extract the information from the background and to profit 
of highly refined interferometric techniques we would have to modulate the 
signal. This result could be achieved, for instance, steadily rotating the whole 
interferometer in the horizontal plane. This solution would however intro- 
duce a new source of noise too, in form of vibrations. It would be better to 
think of a static interferometer with a couple of rotating beams. In practice 
one could have two cylindric, coaxial mirrors; the internal one should be par- 
tially transparent. On the axis one would have a compact rotating head with 
a couple of source/receivers shooting two light beams orthogonally to each 
other; this result would actually be performed with an appropriate beam 
splitter on the axis and a primary source sending light along the axis. The 
couple of cylindrical refiecting surfaces would act as a Fabry-Perot device; 
the return beams would interfere on the axis. 

Of course in this configuration one would have the sought for signal, mod- 
ulated at a frequency double of the rotation frequency of the beams. Any 
imperfection of the mirrors and of the whole set up would of course generate 
perturbations at the same fundamental frequency as the one of the signal. 
The difference between signal and rotation induced noise would be that the 
signal would be peaked on the East- West direction, whereas the various dis- 
turbances would have random orientations of their axes. Repeating many 
measurement runs with different angular configurations of the interferome- 
ter and carefully analyzing the data should allow separating the East- West 
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peaked signal from the rest. 

Besides the rotation induced perturbations one would expect also elastic 
vibration and thermal noise in the solid body of the interferometer. The 
elastic vibrations can be controlled carefully designing the structure in order 
to have proper frequencies not coinciding with the rotation frequency of the 
beams. For a rigid configuration at the scale of the meter, proper frequencies 
are easily greater than 1 kH. 

As for thermal noise, it can be controlled operating at low temperature. 
In any case if a phase difference can be achieved in the order of one hundredth 
of a fringe or better, the signal would easily be bigger than the amplitude of 
the thermal noise of an even moderately cooled device. 

In order to cope with the needs both of mechanical and thermal stability, 
and considering that we need also to properly guide the optical beams a 
best suited material could be sapphire, whose properties make it extremely 
interesting for interferometry in gravitational waves detection |jT2[ . 



The one described above is a simple scheme, showing the principle feasi- 
bility of an experiment. A careful analysis of the technical details would of 
course be needed in order to proceed further. 

We can conclude that the calculations we have written in this paper fix 
the order of magnitude of effects depending from the a? of the Earth and 
show that they should be big enough to be measurable by interferometric 
techniques. 
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